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Q \ Abstract: In this paper we will study quantum field theory of fluctuation modes 

around the rolling tachyon solution on non-BPS D-brane effective action. The goal of 
this paper is to study particle production during the decay of non-BPS D-brane and 
D \ explore possible relation with minisuperspace calculation. We find that the number of 

particles produced on half S-brane exponentially grows for large time which suggests 
that linearised approximation breaks down and also that backreaction of fluctuation 
^ ' field on classical solution should be taken into account. 
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1. Introduction 



A spacelike brane, or S-brane is almost the same as ordinary D-brane except that 
one of its transverse dimensions includes time. S-brane can be also seen as time- 
dependent, soliton-like configurations in a variety of field theories. In string theory, 
the potential for the open string tachyon field on the world-volume of unstable D- 
brane leads to S-branes in a time-dependent version of the construction of D- 
branes as solitons of the open string tachyon. These S-branes can be thought as the 
creation and subsequent decay of an unstable brane. This process has recently at- 
tracted much attention 0, |, |, |, |, |, ^, |, |10|, |12|, |l|, |2|, H, ||, |T], |7l ||, |Tl|, ||, |0|, 

I ^. S-branes have been also extensively stud- 



31, 32, 33, 34, 35 



m 42 



ied in supergravity approach with potentially interesting cosmological applications 

[@, n, n, 13, n, n, g, |i], 0, n, g, m, 0, 0, n, n, 0, |i], ^ n, 0, n, 0]. 

Very nice world-sheet construction of S-branes was given in the classical Qs = 
limit by A. Sen 0, ^, Q where he introduced class of models in bosonic string 
theory obtained by perturbing the fiat space c = 26 CFT with the exactly marginal 
deformation 



S, 



baud 



A / dr cosh X^(t) 



XI] 



where X° is time coordinate, t is a coordinate on the world-sheet boundary and A is a 
free parameter in the range > A < ^ This is family of exact solutions of classical 



^For the most recent discussion of S-branes and other time-dependent processes in string theory, 
see H |8|, H, H, [ttI |78[ |o[ ||. 



^We work in miits a' = 1. 
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open string theory whose space-time interpretation is that of an unstable brane being 
created at a time ~ — r and decaying at a time ~ r with r = — log(sin('7rA)). 

As was shown in [|, ^ this time-dependent process of tachyon condensation 
has many intriguing properties. In particular, it is known that in time-dependent 
backgrounds there is in general no preferred vacuum and particle production is un- 
avoidable. In the open string vacua on S-brane were studied and it was shown 



that they have very mysterious properties. In particular, for ( |1.1|) there is open 
string pair production with a strength characterized by the Hagedorn temperature 
Tr = ^ B 0- As was argued |, ^ this temperature arises from the periodicity 
of the boundary interaction ( |1 . 1|) in imaginary time. 

Recent intensive work on the dynamics of unstable D-branes in string theory 
has led to an effective action for the open string tachyon T and massless open string 
modes (the gauge field on the D-brane) and (the scalar fields describing the 
location of the D-brane in the transverse directions) p|, ^ |^, jl^, |10|, |ll|, |12|, |14, 

As it was shown in 0, ^ the effective action for non-BPS 
21) reproduces several non-trivial aspects of open string dynamics. 



15, 16, 17, 18, 19 



D-brane ([371 



for very nice recent discussion, see [|], |T2|, |T^, |T^, |Tl|, [T5[ . These successes lead one to 
believe that the action ( |3.1| ) captures some class of phenomena in the full classical 
open string theory. As was argued in |jl2[ this action should be thought of as a 
generalization of the DBI action describing the gauge field and scalars on the 
brane ^. On the other hand we should be more careful when we speak about the 
notion " low energy effective action" in case of unstable D-brane. The fact that the 
tachyon has negative m? of the order of the string scale, the notion of an effective 
action which is normally considered as a result of integration out the heavy modes for 
describing the dynamics of light modes, is somewhat unclear here. As was recently 
argued in |]T5| , the situation is even worse by the fact that there are no physical states 
around the tachyon vacuum and hence the usual method of deriving an effective 
action-by comparing the S-matrix elements computed from string theory with those 
computed from the effective action-does not work. The question of usefulness of the 
effective action for non-BPS D-brane was extensively discussed in where it was 
argued that its usefulness can also be judged by comparing the classical solution of 
the equation of motion derived from the effective action with the classical solutions 
in open string theory which are described by boundary conformal field theory. In 
this respect the effective action constructed in [|, has had some remarkable 

success. 

Motivated by this success of effective action description of the open string phe- 
nomena on unstable D-brane and also very nice analysis performed in |]^, ^, we 
asked the question whether there should be some connection between effective action 
description of fluctuation fields around S-brane, or half S-brane, and minisuperspace 



^For review of DBI action, see 
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description of open strings on the same background. In other words, the goal of 
this paper is to construct quadratic action for fluctuation field around half S-brane 
solution on non-BPS D-brane. Since the classical solution describing half S-brane 
explicitly depends on time we can expect that metric and mass term in the action 
for fluctuation field will be function of time as well. As is well known from the study 
of quantum field theory in curved space-time ^ time-dependent background generally 
leads to the creation of particles during the time evolution. Then we could expect 
that particles will be produced during unstable D-brane decay as well. In order to 



study this process we will follow |7^, |7^, |63, |68 |. 

Le us sketch the main idea of particle creation on half S-brane. We start with non- 
BPS D-brane in unstable minimum of the tachyon potential in the asymptotic past 
t —>■ — oo. This corresponds to the situation when the vacuum value of the tachyon 
is zero and small fluctuation modes oscillate around it. Half S-brane solution then 
describes process when the homogeneous tachyon field starts to roll from unstable 
minimum of the potential (T = 0) at the time t = —oo to stable one (T = oo) at 
asymptotic future t = oo. This rolling tachyon solution forms background for the 
fluctuation modes which is manifest-ally time-dependent. It is important to know 
whether these fluctuations, that were small at the beginning of the rolling tachyon 
process, remain small during time evolution or if they become large and hence they 
could affect the classical evolution. Since large value of fluctuation field can be in- 
terpreted as creation of large number of particles we will be mainly interested in 
the time evolution of this quantity. We find that number of particles created during 
D-brane decay exponentially grows at asymptotic future. The similar behavior of the 



fluctuation modes was observed in |3^, ^ , where it was argued that the growing 



of fluctuation modes quickly destabilize the linearised analysis and hence backreac- 
tion of fluctuation field on homogeneous solution should be taken into account. The 
result that we get in this paper also seems to indicate that the linearised analysis of 
fluctuation modes is only suitable for the study of fluctuation at the beginning of the 
rolling tachyon process. This is consistent with the fact that only at the asymptotic 
past where the tachyon is in its unstable minimum open string perturbative states 
exist and one can define open string S-matrix. It is also important to stress that the 
analysis given in this paper is restricted to the pure classical approximation where 
we consider the string coupling Qs equal to zero. However it was shown recently in 
^ [77| , |75| that the coupling between closed and open strings in the rolling 
tachyon process plays fundamental role which suggests that time-dependent D-brane 
decay is very complex problem and its study could answer some fundamental ques- 
tions considering nature of the string theory. 

This paper is organized as follows. In the next section (|^) we briefly review 
minisuperspace analysis of half S-brane to show that particle production during D- 



spor review, see pg|, ITOI, ItI 
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brane decay is natural process from the point of view of quantum field theory on 
D-brane. In section (j^) we introduce a non-BPS Dp-brane effective action proposed 
in ||21|, 0. Section (^) is devoted to study of fluctuation modes above half S-brane 
solution on non-BPS D-brane effective action. In conclusion we outline the results 
obtained in this paper and give their possible interpretation. 



2. Review of the minisuperspace approximation 

In this section we present a short review of minisuperspace approach to the study of 
S-brane dynamics, following seminal papers [H H, 0- 

We wish to understand the dynamics of the open string world-sheet theory with 
a time- dependent tachyon 

S = -— I d^ad'^X^'daX. + [ dr m^iX") . (2.1) 

47r JSa JdT.2 

For the open bosonic string = T where T is the space-time tachyon, while for 
the open superstring ~ after integrating out world-sheet fermions. We use 
the symbol to denote the interaction because the coupling (among other effects) 
impacts a mass to the open string states. We consider three interesting classes 
described by the marginal interactions 

mliX') = ^e^' (2.2) 

m^{X') = ^e-^' (2.3) 
m|(X°) = AcoshX° . (2.4) 

The first case describes the process of brane decay, in which an unstable brane 
decays via tachyon condensation. The second case describes the time-reverse process 
of brane creation, in which an unstable brane emerges from the vacuum. The final 
case describes an S-brane, which is the process of brane creation followed by brane 
decay. Brane decay can be thought of as the future (past) half of an S-brane, i.e. 
as the limiting case where the middle of the S-brane is pushed into the infinite past 
(future). 

In this problem has been studied using the minisuperspace analysis in 

which the effect of the interaction is simply to give a time-dependent shift to the 
masses of all the open string states. In the minisuperspace approximation only the 
zero-mode dependence of the interaction m^(X°) is considered. In this case we can 
plug in the usual mode solution for the free open string with oscillator number N to 
get an effective action for the zero modes 

S = (^-^i;% + (X-l) + 2m2(x°)^ . (2.5) 
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This is the action of a point particle with a time-dependent mass. Here x'^{t) is the 
zero mode part of X'^{a, r), and the second term in is an effective contribution from 
the oscillators, including the usual normal ordering constant. From upper action we 
can write down the Klein-Gordon equation for the open string wave function x), 

[d^d^ - 2m\t) - (iV - 1)) x) = 0, (2.6) 

where [t, x) are the space-time coordinates corresponding to the world-sheet fields 
(X°, X*). This is the equation of motion for a scalar field with time-dependent mass. 
Quantum field theory of the scalar field with time- dependent mass is similar with 



with the QFT in time- dependent background |6^, |7y, |71|] hence we we should make 
few remarks about such field theories. Time translation invariance has been broken, 
so energy is not conserved and there is no preferred set of positive frequency modes. 
As a consequence this leads in generally to particle creation. The probability current 
Jn = ^(0*<9/i0 — <9^0*0) is still conserved, allowing us to define the Klein-Gordon inner 
product 

{f\g)=z jji:^{rd,g-d,rg) , (2.7) 

where S is a spacelike slice. This norm does not depend on the choice of S if / 
and g solve the wave equation. Normalized positive frequency modes are chosen 
to have (/|/) = 1. Negative frequency modes are complex conjugates of positive 
frequency modes, with (/*|/*) = —1. There is a set raising and lowering operators 
associated to each choice of mode decomposition - these operators obey the usual 
oscillator algebra if the corresponding modes are normalized with respect to ( p.7| ). 
We also define a vacuum state associated to each mode decomposition - it is the 
state annihilated by the corresponding lowering operators. 

To illustrate this idea we give an example of the half S-brane corresponding to 
the decay of unstable D-brane so that (|2.6| ) contains following time-dependent mass 
term 

<{t) = ^e* . (2.8) 

Since in the past t ^ — oo we have ordinary D-brane with perturbative open string 
spectrum it is natural to consider \in) vacuum as a vacuum with no particle present. 
Expanding field in plane waves 

0(t,x) = e^P''M(t) (2.9) 

the wave equation becomes 

{d'^ + Xe^ + uj^)u = ,uj^ =p'^ + N -1 (2.10) 
This is a form of Bessel's equation. It has normalized, positive frequency solutions 



A-^^^^^i^J_2..(2v^e*/^) (2.11) 
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where superscript in and out on a wave function denotes solutions that are purely 
positive frequency when t— ooort— i>oo. Using properties of Bessel's functions 
we can easily show that m*" approaches flat positive frequency plane waves in the far 
past if: — > — oo 

M^" ~ -Le"*"* . (2.12) 
We can also consider the wave functions 



^out 



y|(,e-^)-i/2^(2)^(2v^,*/^) , (2.13) 



that are purely positive frequency in the far future t ^ oo 

^out _ _ 2^^/Xe'/^} . (2.14) 

V z 

Generally and u'^^ are related by celebrated Bogolubov transformations 

^out ^ j^^in ^ ^^.n* ^ ^2.15) 

where coefficients A,B can be determined from known solutions In partic- 
ular, for modes given above it can be shown that 

T • (2.16) 

sinh27rcjr(l - 2icu)y ^ ' 



These coefficients obey unitarity relation — = 1. The relation ( |2.15| ) between 
in and out modes implies the relation between in and out creation and annihilation 
operators 

ttin = Aaout + 5*a),„i . (2.17) 

From this relation we see that the condition aj„ \in) = implies that \in) is squeezed 
state 0. Physically, this is the statement that particles are produced during brane 
decay: if we start in a state with no particles at t ^ — oo, there will be many particles 
at time t —* oo. The density of particles of momentum p is given 

np = |7|2 = e"^"'" . (2.18) 

We see that despite the fact that \in) is pure state this is precisely the Boltzmann 
density of states at temperature = l/47r. In string units, is Hagedorn tem- 
perature. As was argued in |^, |^ the fact that this temperature is so high means 
that Qs corrections are likely important even for Qs ^ 0. 

The main conclusion which we can get from this brief review of minisuperspace 
analysis of S-branes is that in the limit Qs ^ there is a particle production during 
D-brane decay. We would like to see the similar behavior from the effective action 
point of view. We will discuss this problem in next sections. 
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3. Non-BPS D-brane effective action 



The dynamics of the tachyon field T on non-BPS D-brane can be described by 
effective field theory where the value of T and its derivatives satisfy some conditions 



ri|, |T2| |13|. The effective action for real tachyon T on non-BPS Dp-brane in Type 



II string theory is expected in the form |T^, |18|, |I9| 

S = J (F+^xL = -MpJ d^^^xV{T)^l + ^r%TdJf , (3.1) 
where Mp is tension of non-BPS Dp-brane ^. Another form of the effective action for 



non-BPS D-brane has been obtained recently in ||12 



S = j (F+^xL = -Mp J (F+^xV\T)Jl + ^ + r]t^-d^Td,T , V^{T) = ^ 



^ 2 



(3.2) 

where the tachyon field T given in ( p. 21 ) is related to T in ( |3.1| ) by the field redefinition 



T r 

—pz = sinh ^= . (3.3) 
V2 V2 ^ ^ 



We will concentrate in this paper on the action ( |3.2| ) since it has the tachyon solution 
of equation of motion which has the same form as classical solution in open string 
theory that in turn is described by boundary conformal field theory (BCFT). For 
that reason we believe that field theory for fiuctuation modes in the effective action 
( |3.2|) could have some relation to the quantum field theory reviewed in the previous 
section that arises as minisuperspace approximation of exact BCFT. 

We would like also study the behavior of massless modes on non-BPS D-brane 
effective action during D-brane decay. For that reason we should include them in 
In order to do this we use an equivalence between ( |3.2D and ( ^.1[ ) where the 



action (p.l|) containing massless modes is well known [|l^, jl^, [H, [T9| 



S= -Mp 



J (F+\V{T)^- det(r/^, + d^Y^d,YJ5u + F^, + d^Td^T) , (3.4) 



where and for yU , z/ = 0, . . . , p , / = p+l, . . . , 9 are the gauge and the transverse 
scalar fields on the world-volume of the non-BPS D-brane. If we expect that the field 
redefinition (|3.3|) remains the same even with nonzero massless fields and also that 
these massless modes are the same in both forms of the action we can write 

T = y2sinh^ ,v(r) = — ^ = — ^ — , d.r = -ML_ (3.5) 



are using the convention where fj^i, = diag{-l, 1, . . . , 1) , /i, = 0, . . . ,p , a, 5 = 1, 
and the fundamental string tension has been set equal to (27r)^^ ,«' = !. 
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and hence we obtain generalized form of (|3?2 



where 



S = -MpJ dP+^xV{T)V-detA , 



(3.6) 



A 



V{T) = , F,, = d^A, - d,A. 



1 + ^ 



(3.7) 



Action (|3.6|) will be useful when we will discuss the massless mode on half S-brane. 
For classical solutions that depend on T only the action (|3.2|) is sufficient so that the 
equation of motion that arises from it is 



T 



T 



2(l + f)VB 



'1 



T2- 



where 



B = 1 + ^ + v'^'d^Td^T . 



Solution of ( |3.8| ) that is interpreted as half S-brane is given by 



(3.8) 



(3.9) 



(3.10) 



where a is some constant which can be set to 1 by time translation! 21 



In the next section we present an analysis of the fluctuation modes above classical 
solution ( |3.1CI| ). 



4. Analysis of fluctuations 



In this section we will analyse the fluctuation field around the classical half S-brane 
solution (|3.1CI|) . To do this we write the tachyon field T that appears in ( |3.2|) 



as 



T(t,x) = r,(t) + 0(t,x) 



(4.1) 



where Tc is given in ( |3.1U| ) and (f){t, x) is fiuctuation field around where we presume 
that it is small with respect to the classical solution. 

r2D and perform an expansion with respect to 



We insert ( |4.1| ) into ( p.2D and perform an expansion with respect to (f) up to 
second order. As a result we obtain a quadratic action for free massive scalar field 
(j) where metric and mass term are functions of classical solution Tc. More precisely, 
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let us write 



S = - J dxL{Tc + (/ 
(6L{T, 



— dx 



-- I dx 



bTbT 



+ 2 



6T5d^T' 



+ 



\ 5T 
6d^T5d,T 



+ 



dxL{Tc 



d„ 



+ ... 



(4.2) 



where dots mean terms of higher order in 0. Using equation of motion and integration 
by parts we can easily show that the expression on the second line in (|4^ ) is equal 
to zero. The quadratic effective action for is then given on the third line in 
We can rewrite it in more familiar form using 



dx 



5T5d^T' 



dxd,. 



5T5d^T' 



dxdu 



dx 



6T6df,T 



dx 



5T5dy 



dxdu 



ST5d^T 



(4.3) 



where as usual we do not carry about boundary terms. Then we obtain the action 
for fluctuation modes up to second order in 



S = — / dtd^K 
2 



G-(Te) = ^Sii,M^(T.) = S-5„ 



Sd^TSd,.T 



ST5T 



6T6df,T 



(4.4) 



For half S-brane solution ( p.lO| ) we obtain from ( |3.2|) following metric and mass term 
that appear in ( [4.4| ) 



G' 



00 



MpV^ MpV'^T^ 



(B)3/2 



a 



ab 



M 1/2 (^"^ M S"-'' 



ab 



1 + 



M= = -Ml 

2 



(4.5) 



where for ( p.lO| ) we have B = l + ^ — = 1. According to ( ^4.5| ) the action for 
fluctuation modes around half S-brane solution has the form 

1 „ , , 1 



S = I dtdyi 

2 



-a 



+ 



:,V2t 



(4.6) 
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This form of the action will be starting point for the analysis of particle production 
on half S-brane. We will show that ( |4.6| ) has similar form as the action describing 
scalar field with time-dependent mass. Such quantum field theories are well known, 
for very nice and detailed discussion, see in particular |6^ . 

To begin with, we should determine the conjugate momentum to (j) and the 
Hamiltonian that arise from (14. 61) 



n(t,x) 



5L 



n{t,:>i) = n(t,x)0(t,x) - L 



1 

2 [JL 



+ G'"'(t)da(pdb(p + M^'J 
Hit) = J dxH{t,^) 



(4.7) 



Since the effective metric and mass are functions of time only it is natural to write 



operator 0(t, x) as 



0(t,x) 



dkQuit)e^^'' ,n(t,x) = Mp I rfkPk(t)e 



ikx 



(4. 



where we have included the normalization factor 7;^- into the definition of dk. When 
we insert ( [4.8|) into the Hamiltonian we get 

M 



j dk[PkP-k 



G{t)k^ + 



M2 



(4.9) 



Condition (pit, x) = 0^(t, x) implies = Q-k and consequently the Hamiltonian for 
the scalar field can be seen as collection of Hamiltonians i^k of quantum oscillators 
with complex coordinates Qk 



H = MpJ dkH^ ^H^=l {P^Pl + ^lit)Q^^Qi 



and with the time- dependent frequencies 

/Vf2 
IVlp 

Canonical equal time commutation relations for scalar field 

[0(t,x),H(t,y)] = z5(x-y) , 
[0(t,x),0(t,y)] = , 
[H(t,x),H(t,y)] = , 



6 • 



(4.10) 



(4.11) 



(4.12) 
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imply following commutation relations between Qk, -Pk 

[Qk, Pk'] = ^S{k + k') , [Qk, Qk'] = [Pk, Pk'] = . (4.13) 

In the following we will review the basic properties of quantum oscillators with time- 
dependent frequencies. We will mainly follow . 

We would like to go to a time-dependent frame in Hilbert space in which the 
Hamiltonian ifk is diagonal at every moment of time. For that reason we define 
time-dependent operators ak{t) and al_{t) by 

Ok = (Qk^k + ^^k) , 

(4.14) 

These operators are Hermitian conjugate and have the standard commutation rela- 
tion of creation-annihilation operators 

[a^,al] = 6{k-k') . (4.15) 

Using these operators one can write 



y^L^^Jn^dt^^_^iJn^dt^Q 



(4.16) 



so that the Hamiltonian i^k has the form 



i^k = ^ (l + 4«k + a^-k«-k) • (4.17) 

The orthonormal frame in Hilbert space which diagonalizes the Hamiltonian ifk is 
given by acting creation operators aj^ on time dependent vacuum state |0j) which is 
annihilated by the operator ak(t). The operators ak(t) , al_{t) obey following equation 
of motion 

A, - ilLp2i / n^dt \ 

- 2^]k 

. + ^k 



''^-2nk 



(4.18) 
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The solution of ( |4.18| ) can be written in terms of constant creation and annihilation 
operators Oq ^ , ao,k as follows 



aLi,(t) = /9k(t)ao,k + al(t)a, 



t 

0,-k ■ 



If we insert (|4.19| ) to ( [4.18| ) we easily get 



^k - 

e 



2n 



(4.19) 



(4.20) 



System ([4.1 9|) represents so named Bogolubov transformation between two pairs of 
creation-annihilation operators. If the oscillator initially {t = — oo) is in the vacuum 
state then its state |Oo) is annihilated by the operator ag.k and the initial condition 
for functions a^it) , Pk(t) are 



|ak(-oo)| = 1 ,/9k(-oo) = . 



(4.21) 



At the moment t it will not be in the vacuum state |0t) annihilated by the operator 
ak(t). Recall that we are working in Heisenberg representation where the state does 
not evolve. There is relation between the states considered 



ioo) = n 



exp 



^|ak(t)| V2«kW 



(4.22) 



At the moment t the average number of particles produced in the k-th mode is 



A^k(t) = (Oo|aL(t)ak(t)|Oo) = |/3k(t)|' • 
In Heisenberg representation the equations of motion for Qk , -Pk are 

Qk = t[H^, Qk] = a , a = -^l{t)Qk . 

Their solutions can be written using operators ao,k , clq _k as follow 
where modes ,Q^^^*{t) = Q^^^ obey differential equation 

Qt^ + nimi^^ = . 

We can easily express /3k (t) in terms of these modes 



/9k(t) 



V2n 



(4.23) 



(4.24) 



(4.25) 



(4.26) 



(4.27) 
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so that the number A^k(^) of particle produced is equal to 




(4.28) 



Now we apply the formalism reviewed above to the case of fluctuation field around 
half S-brane solution ([4.61). Since the number iVk(t) of particle produced ( [4. 28] ) is 
function of obeying ([4.26|) , we must solve this equation. In case when we do not 
need to know an exact solution the most efficient method to solve ([4.26|) is to use 
WKB approximation (For review of the application of this method in QFT in curved 
space-time, see for recent discussion, see [0.) In fact, we are mainly interested 
in estimation of A^k(^) at asymptotic future so that the WKB approximation can be 
applied. Setting 



Q 



(-) 



1 



:exp(-z J dtW]i) 



one can see that the solution of equation ( [4.26|) is equivalent to the condition 



If following condition holds 



< 1 ,Qk 



SWT 




(4.29) 



(4.30) 



(4.31) 



we can write 



W^{t) = a(t) . (4.32) 

Let us apply WKB approximation for equation ( |4.26| ), following very nice analysis 
given in [|73|. Firstly we observe that there is a time t^, when vanishes 

2 _ 1.2 ^ 



ni{Q = ^ e 



(4.33) 



i.e. t=K is classical turning point. According to ||7^ we define region I as the region 
such that > and region II the region where fi^ < 0. At some time t deep in 
region I the WKB approximation of can be written as 

A 



Q 



(-) 

k 



expi^— « 



dtQh 



(4.34) 



\/2f2k -^ti 

where ti is the initial time at which the normalization is performed. We can also 
easily see that deep in region I the condition ( [4.31| ) is valid. Since for t — — oo 
the half S-brane solution (|3.10|) tends to zero corresponding to the original non-BPS 
D-brane it is natural to demand that for t — oo approaches the plane wave 
mode 

fr. exp(-?c^o,k(i 

V ^^0,k 



ti)) 



(4.35) 
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In fact, deep in the region I we have 



I 1 z 

^o,k - k e 



(4.36) 



and hence ■* is equal to 



A -«a;o,k(t-ti)+»77| 



(eV^t-eV^ti) 



(4.37) 



Comparing with (|4.34| ) we immediately obtain the value of normalization constant 
A equal to 1. We must also stress that we are considered the scalar modes with the 
initial frequency cUq = — ^ > since these correspond to the standard fluctuation 
modes on non-BPS D-brane in the past infinity. For modes that obey — ^ < the 



analysis is much more complicated, as was shown in |75, 76 



In the region II the situation is different. In this case f2k becomes complex. 
However this does not prevent to use WKB approximation |]74|. If we write nk(t) as 
^k(^) = ^|^k(i^)| then the solution in region II in WKB approximation is 



|a(t)|i/2 



exp 



rft'|fik(t')| 



+ 



ia(t)r/' 



exp 



rft'|fik(t')l 



(4.38) 



This solution should be matched on the solution in the region I. At the turning 
point t = the WKB approximation breaks down and we must use the usual WKB 
procedure and approximate the potential around the point t^, by a straight line 
such that 

nl{t) ^ -a{t - U) , (4.39) 



where 



a 



dt 



1 + 2u;l^ 



> . 



(4.40) 



The solutions of the equation of motion with such a potential are given in terms of 
Airy functions of first and second kinds 



5iAi(s) + B^Biis) , s = a^'\t - Q . 



(4.41) 



Now we use the asymptotic behavior of the Airy functions to calculate the relation 
between Bi , B2 with A on one hand and C± on the other. In region I, for a value of 
t not too far from t*, ( ^.41| ) can be written as 



a 12 



{Bi - iB2) exp 



2i 
— ( 
3 

2i 

r 







3 




\t 




2 


+ t) 












3 




\t 




2 











+ 



(4.42) 
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while in region II, under the same conditions, the function Q^. ^ (t) can be expressed 
as 



_ j_ 

a 12 
2V2 



2 1 , ,3 



(4.43) 

Now we evaluate in region I the solution ([4 .341) for the potential ([4.39|) . The integral 
in the exponent can be written as 



(4.44) 



The frequency ( |4.39| ) is used in the second integral only, assuming that t is not too 
far away from t^,. The quantity $k is just a number and its calculation would require 
the knowledge of f2k in the whole region I. However as we will see it does not enter 
in final result and therefore we are not interested in its value. Then we find 



Qk,i (^) = « ^1^-^*1 4exp ( ya2|t-t|2 
so that comparing with ( [4.42| ) we find 

Bi = iB2 ,-82 = v^a^^e*''^? . 



(4.45) 



(4.46) 



In the same way we can establish the link between C+, C_ and Bi,B2. This can be 
done in the same way as in the region I. If we define \l/k as 



\^k{t')\dt' 



(4.47) 



then ( [4.3^ ) reduces to the expression 

Q^-j)(t) = a-^t-Q--*\ C+e^^ exp 



exp 



--«5(t _t^)5 



+ 



By comparison with ( [4.43| ) we get 



2v^ 



a 



TT 



(4.48) 



(4.49) 



The upshot of this analysis is that generally both coefficients C+, C_ are nonzero 
when we match the solution in region II to the solution in region I that corresponds 
to the plane wave at asymptotic past t ^ — 00. For t ^ U the leading term in qI_ "* 
is proportional to C_ so that we have 

Bo 



v/?|2a 



- JJ dt'\nut')\ 



(4.50) 
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and consequently the number N]^{t) of particle produced for large t is 

Ar^(t) = ^e-^/.*''^*'l^''(*')l--. (4.51) 

In asymptotic future t tj and we get 

N^{t) ~ e'*'' ~ e^'f , (4.52) 

where we have used the fact that for t ^ oo This result implies that 

during D-brane decay infinite number of particles is produced. The result is in agree- 
ment with p8| , p9| , , where similar exponential instability of fluctuation modes was 
found in case of different non-BPS effective actions. It was then argued there that 
this result implies that the approximation, where we consider fluctuation modes as 
small with respect to classical condensate and hence we can neglect their backreac- 
tion, breaks down very early from the beginning of the time evolution. We observe 

the same behavior in our approach as well. More precisely, we can easily see that 

t 

when we insert asymptotic form of fluctuation ~ e to the action ([4.6|) that this 
action dominates over the action (13.21) evaluated on half S-brane solution 



SiT,) = -Mp J dt^^^ . (4.53) 

Then it follows that the approximation when we regard as small fluctuation modes 
without any backreaction on classical solution breaks down at late times. It would 
be nice to perform analysis of the fluctuation modes around half S-brane solution 
when we take this backreaction into account. We hope to return to this problem in 
future. 

4.1 Massless fluctuation on half S-brane 

In this section we extend our analysis to the case of the massless fluctuations on half 
S-brane. We start with the action 

S = -MpJ dP+'xV{T) sJ- det (r/^, + d^YW,Y-^5ij + F^, + V\T)d^Td,T) . 

(4.54) 

We will focus on scalar modes only since for gauge fleld the situation is 
completely the same. When we perform an expansion of ( [4.54|) around classical 
solution Tc up to second order in Y^ we get 

S = -Mpj dP+'xV{T,)^- det {r]^, + V^{T,)d^T,d,T, + d^Y^d^Y^du) = 

= -MpJ d-P+^xViT,)^- det g{T,)^ det^d'i + g^^^^d^Y^d^YJ 5ij) = 

= -MpJ dP+'xV{T,)^~detg{T,)- 

/ dP^'xViT,)^-detgin)g^''iT,)d,Y'd,Y'5jj , 

(4.55) 
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where 

Q{Tc)^.u = V,. + V\T,)d^TAT, (4.56) 
which for half S-brane solution ( p.lC ) is equal to 

g{T,)^, = diag I YTW^ , 1, • • • , 1 I • (4.57) 

Consequently the quadratic effective action for scalar fluctuation modes around 
half S-brane solution has the form 

S = / d^+'xVin)^~ det giT,)g>^''in)d,Y'd,Y'6jj = 

= -^j d^dt [-d,Y^d,Y, + ^^^d.Y^d,Y,j . 

(4.58) 

Let us consider one scalar mode, say = Y. From ([4.58|) we get conjugate 

momentum to Y 

U{t, x) = -J^ = M,Y{t, x) (4.59) 
dY[t, x) 

so that the Hamiltonian is 



H = ^ J d^{UY - l) = ^ J d^(^^ + M,G{tr'daYd,Y 



(4.60) 



If we take take y = / (iklk(t)e*^^ the Hamiltonian reduces to the collection of 
Hamiltonians ifk with frequency ^\{t) 

H = JdkH^,H^ = \ + nl{t)Y^Y^^ , (4.61) 

where 

As we know from the previous section in order to determine number A^k(^) of particle 
produced at time t we should solve following equation 

Qi~^ + nl{t)Qi-^ = , (4.63) 

where f2k is given in ( [4.621) . As opposite to the case of tachyon fluctuation mode 
studied in previous section from fik given in ([4.62|) we get 

lim — = — \= , lim nJt) = (4.64) 
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and hence WKB approximation cannot be used. However in order to estimate num- 
ber of particles at asymptotic future we does not need to know an exact solution in 
the whole region, rather we could try to solve (|4.63|) in asymptotic future. For large 
t we can write 

nl{t) = ee-^' (4.65) 

so that (l^eD is 

q(-) + k^e-^'Qi^^ = . (4.66) 
After substitution m = y/2ke ^ we obtain following equation 

m'^-Skl + rn^^ + m^Q^") = ■ (4.67) 
d'^m dm 

This is Bessel's equation of the first kind with the general solution 

Qi:\m) = AMm) + A2Yoim) . (4.68) 



For m — > we have 



and hence 



2 

Tfl 

Jo(m) ~ 1 - — , Yo{m) ~ Inm (4.69) 

Qi-\t) ~ Ai(l - ye--^*) + A^t . (4.70) 

Let us consider for a moment the situation when A2 = 0. Then the number A^k of 
particle produced for large t is equal to 

A^k ~ f^k — n h A? A? , ~ — (4.71) 

\ 2 nl 2 ^ e-73 2 2 ^ ' 

which is clearly nonphysical result since the number of particles cannot be negative. 
For that reason any physical solution must contain term proportional to A2. However 
this term dominates at large time so we can write 

gir^ ~ A2t . (4.72) 

For such a solution the number A^k(^) of massless particles produced during D-brane 
decay is equal to 



A^k ~ f ^ + A^t'j - - ~ ^ ~ . (4.73) 



We see that the number of particles that are created during the D-brane decay 
exponentially grows. This is consequence of the fact that for t ^ 00 the effective 



frequency goes to zero. Following |38, B9, 40 1 we can interpret this exponential 



growth of particles as breaking of linearised approximation and their backreaction 
on classical solution should be taken into account in the process of D-brane decay. 
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5. Conclusion 



In this paper we have studied the particle creation during the time-dependent process 
of non-BPS D-brane decay (half S-brane) in the effective field theory approximation. 
We hoped that the behavior of fluctuation modes around rolling tachyon solution on 
non-BPS D-brane could reflect the behavior of open string modes as was described 
using minisuperspace approach ^, 0]. In particular, we have been mainly inter- 
ested in the question whether particle production occurs during D-brane decay and 
whether vacuum states of fluctuation modes in asymptotic past and future have sim- 
ilar properties as vacuum states introduced in f^, |, . While the particle production 
was unavoidable in our effective action description as well since the action for fluc- 
tuation field around half S-brane solution corresponds to the quantum field theory 
action in time-dependent background, we have find that linearised approximation is 
not completely sufficient to reproduce results given in |0, ||, . In particular we have 
shown that number of fluctuation modes exponentially grows for large time which 
implies that the linearised approximation breaks down. Comparing the effective ac- 
tion for fluctuation field given in this paper and the quantum field theory action given 
in |, 1^ we immediately see reason why our results are so different. In p, ^, 
effective frequency for given mode grows to infinity in the asymptotic future and 
hence given mode exponentially decays. On the other hand, in the quadratic action 
for fluctuation field given in previous sections the time-dependent frequency either 
approach constant value, as in the case of tachyon fluctuation, or approaches zero for 
the case of massless fluctuations. This behavior of effective frequency results into the 
exponential growth of number of particles produced during D-brane decay. Similar 
exponential instability of fluctuations has been observed in |3^, where this 
huge growth of number of particles produced has been interpreted as breaking of the 
linearised approximation. To conclude we mean that the linearised description of 
fluctuation modes around half S-brane is only efficient at the beginning of the rolling 
tachyon process. After some time it breaks down and we should take into account 
backreaction of fluctuation field on the classical solution. 

We must also stress that in this paper, following H, 0, we restrict to the case 



gs = 0. However recent results given in very interesting papers |7^, [79 



show that in order to correctly describe D-brane decay we should take into account 
interaction between closed and open strings and hence to go beyond the free case 
= 0. 

To conclude, we hope that our modest results given in this paper could be helpful 
for other study of non-BPS D-brane effective action and time-dependent process of 
D-brane decay. 
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